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The early stage of high multiplicity nuclear collisions is represented by a nearly quarkless, hot, deconfined pure gluon 
plasma. This new scenario should be characterized by a suppression of high pr photons and dileptons as well as hy reduced 
baryon to meson ratios. We present the numerical results for central Ph-l-Pb collisions at the LHC energies by using the 
ideal Bjorken hydrodynamics with time-dependent quark fugacity. It is shown that about 25% of final total entropy is 
generated during the hydrodynamic evolution of chemically undersaturated quark-gluon plasma. 
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1 Introduction 

The proper understanding of the initial and the early stage 
of ultra-relativistic pp-, pA- and heavy ion AA- collisions 
is a topic of great importance for our understanding of hot 
and dense QCD matter formed in the laboratory and in the 
early universe. One of the central questions is how the initial 
highly nonequilibrium system evolves to a state of partial 
thermodynamic equilibrium at later stages of nuclear colli¬ 
sions. There exist several models which describe the initial 
state in terms of non-equilibrium parton cascades (Wang & 
Gyulassy 1991; Xu & Greiner 2005), minijets (Eskola & 
Kajantie 1997), color glass condensate (McLerran & Venu- 
gopalan 1994), coherent chromofields (Magas et al. 2001; 
Mishustin & Kapusta 2002) etc. It is commonly believed 
that the strong non-equilibrium effects persist only for 
a short time ~ 1 / Qs, where Qs - 1-2 GeV is the so-called 
saturation scale (Gribov et al. 1983), but at later times the 
system reaches a state of a partial thermodynamic equilib¬ 
rium. 

Relatively large gluon-gluon cross sections lead to the 
idea (van Hove & Pokorski 1975) that the gluonic com- 
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ponents of colliding nucleons interact more strongly than 
the quark-antiquark ones. Then the two-step equilibration of 
the quark-gluon plasma (QGP) was proposed (Raha 1990; 
Shuryak 1992; Alam et al. 1994; McLerran & Venugopalan 
1994; Krasnitz & Venugopalan 2001). In this approach the 
gluon thermalization takes place at the proper time Tg < 
Ifmjc and the (anti)quark equilibration occurs at t* > Tg. 

The surprising result, that only very few soft quarks are 
present at an early stage of a relativistic collision, was ob¬ 
tained in many transport calculations (Biro et al. 1993; Roy 
et al. 1997; Elliott & Rischke 2000; Blaizot et al. 2013; 
Uphoff et al. 2015). Observable consequences of the above 
two-step scenario was considered by several authors, see 
e.g. (Strickland 1994, Kampfer & Pavlenko 1994, Traxler 
& Thoma 1996; Dutta et al. 2002; Gelis et al. 2004; Scar- 
dina et al. 2013; Liu et al. 2014; Monnai 2014; Stoecker et 
al 2015; Vovchenko et al. 2015a). One immediate prediction 
of such a ’’pure glue” initial scenario is suppressed yields of 
hard ’’thermal” photons and dileptonsQ. Such a suppression 
occurs due to the reduction of the electric charge density 


* One should distinguish such particles from photons and Drell-Yan 
dileptons produced in inelastic collisions of initial partons. 
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Fig. 1 (Color online) The quark fugacity /I as a func¬ 
tion of proper time t. The solid and dashed lines corre¬ 
spond to chemically undersaturated (US) plasma with pa¬ 
rameters Tt = 1 fm/c and 5 fm/c, respectively. The dashed- 
dotted line corresponds to the case of chemical equilibrium 
(T = 1). 

as compared to chemically equilibrated quark-gluon plasma 
(QGP) at early stages of the reaction. 

2 Evolution of undersaturated QGP in 
nuclear collision 

Below we assume that a thermally (but not necessary chemi¬ 
cally) equilibrated QGP is created initially in a nuclear colli¬ 
sion. In this section we use the parameters typical for central 
Pb-i-Pb collisions at the LHC energy -yT/v/v = 2.76 TeV. By 
using the one-dimensional scaling hydrodynamics (Bjorken 
1983) we consider the space-time evolution of QGP pro¬ 
duced at the proper time t - tq. We adopt the equation of 
state of an ideal gas of massless gluons, quarks and anti¬ 
quarks. At zero net baryon density one may describe devia¬ 
tions from chemical equilibrium for quarks and antiquarksQ 
by introducing the quark fugacity A. Within these approxi¬ 
mations, the following relations for the energy density s and 
pressure P can be written (h-c-l) 

s^3P^aT\ ^(l6 + /ly/V/j, (1) 

where T is temperature, Nf is the number of quark flavours 
(unless stated otherwise, we assume that Nf - 3). The 
first and second terms in the last equality describe, respec¬ 
tively, the contributions of gluons and quark-antiquark pairs 
to the energy density. The parameter A changes from zero 
for the pure gluonic system to unity for chemically equili¬ 
brated QGP. 


By using the relation P — e/3 and neglecting the vis¬ 
cosity effects one can easily get the analytic solution of the 
hydrodynamic equation in the Bjorken model 



^ In accordance with the two-step approximation (see Sec. 1) we ne¬ 
glect deviations from chemical equilibrium for gluons, assuming that their 
fugacity is equal to unity during the whole process of the QGP evolution 
at r ^ To. 


where the parameter tq corresponds to the initial proper 
time of the hydrodynamic expansion. In principle one may 
determine A and T as functions of t by solving numerically 
the additional rate equation for quark density evolution (see 
e.g. Biro et al. 1993; Monnai 2015). The qualitative analysis 
can be performed by introducing the analytic parametriza- 
tion (Vovchenko et al. 2015a) 

' To - 


/1(t) = 1 - exp 


(3) 


where t, is the model parameter characterizing the quark 
chemical equilibration time. Calculations of different au¬ 
thors lead to different estimates for t, , ranging from 
T* ~ 1 fm/c (Ruggieri et al. 2015) to t, ~ 5 fm/c (Xu & 
Greiner 2005). One should have in mind that this parameter 
may depend on the combination of nuclei and the bombard¬ 
ing energy. We expect that t* will be larger for peripheral 
events and lighter combinations of nuclei. Figure [1] shows 
the A{t) values for different choices of the parameter t, . 

Introducing the quark chemical potential p = T In T and 
using thermodynamic relations, one can write down the fol¬ 
lowing expression for entropy density of the QGP 

'TO 2 

i -HT(0.66-0.161nT)V/]. (4) 

By using ([T]i-(l4ll one can show that st is increasing function 
of T, i.e. ST ^ iflTo, where the equality holds only in the 
equilibrium limit A - 1. 

Within the Bjorken model the total entropy per unit 
space-time rapidity n = tanh“*(z/f) can be expressed as 
(Satarov et al. 2007)0 
dS{T) 


drj 


= nRls{T)T, 


(5) 


where Ra is the geometrical radius of the colliding nuclei. 
Therefore, we obtain that the total entropy per unit space- 
time rapidity is not conserved: it gradually increases during 
the system expansion from the pure glue initial state. Note 
that this increase occurs within the ideal hydrodynamics, 
in absence of viscosity effects. We think that future models 
for extracting the viscosity values from the observed data 
should take into account the suppression of quarks at the 
initial state of a nuclear collision. 

To fix the initial temperature, we assume that the 
Bjorken solution is valid until the freeze-out hypersurface 
T = Tf. The latter is determined by the condition T{Tf) - 
156 MeV. Such a temperature has been extracted (Stachel 
et al. 2014) from the thermal fit of hadron ratios in the con¬ 
sidered reaction. In our calculation we use the approximate 
relation (Hwa & Kajantie 1985) between the total entropy 
per unit space-time rapidity and the rapidity density of pi- 
ons 


dS{Tf) 

drj 


dy 


1.7- 10^ 


( 6 ) 


where v ^ 6.3 is the entropy per pion at freeze-out 
(Vovchenko et al. 2015a) and dNnjdy 2700 is the ob¬ 
served yield of pions at midrapidity (Abbas et al. 2013). 


We consider a purely central collision of equal nuclei. 
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Fig. 2 (Color online) Temperature of QGP produced in 
central PbPb collision ( -sJsnn - 2.76 TeV) as a function 
of T. The solid and dashed curves correspond to chemically 
undersaturated matter assuming the parameters T, = 1 fm/c 
and 5 fm/c, respectively. The dashed-dotted line is calcu¬ 
lated within the equilibrium scenario. 



Fig. 4 (Color online) Mass distribution of thermal dilep¬ 
tons in central PbPb collisions at = 2.76 TeV. The 

solid and dashed lines are calculated with parameters t* = 
1 fm/c, and 5 fm/c, respectively. The dashed-dotted line 
corresponds to the chemical equilibrium scenario. 



Fig. 3 (Color online) Same as Fig.|2] but for the entropy 
per unit space-time rapidity. 


By using this procedure we have calculated the temper¬ 
ature and entropy density of matter in the pure glue initial 
scenario and compared the results with the chemical equi¬ 
librium case. Figure |2] shows the evolution of temperature 
for the same values of parameters t* and tq as in Fig. [1] 
One can see that temperature of the US plasma is notice¬ 
ably higher than in the equilibrium scenario. This increase 
is especially visible at t < t» . 

Figure [3] shows the results for the total entropy evo¬ 
lution. One can see that this quantity gradually increases 
and reaches the freeze-out value (|6]) during the time interval 
At ~ T,. According to our calculations, the total increase 
of entropy is not sensitive to t* and equals about 25% of the 
hnal value. 

3 Dilepton and photon spectra in the pure 
glue initial scenario 

Hard thermal photons and dileptons are sensitive probes of 
hot initial stages of high energy nuclear collisions. We as¬ 
sume that hard dileptons are mostly produced in the decon- 


hned phase via the qq UU annihilation processes. Below 
we again apply the Bjorken hydrodynamics to describe the 
evolution of QGP in a heavy-ion collision. As compared to 
early calculations (Hwa & Kajantie 1985; Kampfer et al. 
1990) we include the additional factor which takes into 
account the reduction of quark and antiquark densities in 
the chemically nonequilibrium case. We get the following 
expression for the mass spectrum of e^e~ pairs: 


dM^ dy 


^ Z 


i~u,d,s 



TdrTij) K\ {MITiT)} A^(t) , 


(7) 


where a - is the electromagnetic coupling constant, q/ is 
the charge of the quark flavour i in units of e and Ki (x) is 
the Macdonald function. As above, current masses of quarks 
are disregarded for all flavours. Note that different scenarios 
considered in Sec. 2 correspond to different choices of /1 (t) 
and T (r). The results of numerical calculations are shown in 
Fig.® One can see that at M > 1 GeV/c^ the dilepton spec¬ 
tra are strongly sensitive to chemical nonequilibrium effects. 
A special investigation shows that such dileptons are mainly 
produced at hot early stages of the reaction. In the pure glue 
initial scenario yields of hard dileptons are suppressed as 
compared to the equilibrium case. 

To study the emission of hard thermal photons we pro¬ 
ceed from analytic formulae for chemically equilibrated 
QGP suggested by Kapusta et al. (1991). In the lowest-order 
approximation in strong coupling constant as main sources 
of the real photon production are the qg and "qg Compton 
scatterings as well as the q'q annihilations. Attempts to con¬ 
sider the chemically-nonequilibrium scenario have been al¬ 
ready made in (Strickland 1994, Kampfer & Pavlenko 1994, 
Traxler & Thoma 1996; Dutta et al. 2002; Gelis et al. 2004). 
Following their procedure, we include additional suppres¬ 
sion factors A and A^ into the components of photon produc¬ 
tion corresponding, respectively, to the Compton scattering 
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Fig. 5 (Color online) Spectrum of the thermal photons as 
a faction of transverse momentum in central PbPb collisions 
at -s/snm = 2.76 TeV. The solid and dashed lines corre¬ 
spond, respectively, to t* = 1 fm/c and 5 fm/c. Dots show 
experimental data (Lohner et al. 2013) for (0^0)% most 
central events. 


and qq annihilation terms. We have obtained the following 
expression for the invariant momentum distribution of hard 
thermal photons (Vovchenko et al. 2015b): 


„ 4q; V 

j2 j “ 3 2 j 

d prdy n i=u,d, 




q^R^ Tdra^T^x ( 8 ) 


Here y and pj- are the rapidity and transverse momentum 
of photons (it is assumed that pr > T), and the constants 
a ^ 0.20, b ^ 0.99, c ^ 0.88. Below we disregard the tem¬ 
perature dependence of strong coupling constant, assuming 
that - 0.3. 

Figure |5] shows the photon spectra calculated for the 
same reaction and same model parameters as above. Again 
one can see a noticeable suppression of high pr photon 
yields as compared to the equilibrium scenarioQ Accord¬ 
ing to Fig. 13 the observed data are better reproduced for 
smaller values of t* . Note, that our calculation does not in¬ 
clude the contribution of photons from initial parton-parton 
collisions. This prediction should be verified in more real¬ 
istic calculations, taking into account the transverse flow of 
deconfined matter. 


4 Evolution of the pure-glue matter 

It is instructive to consider qualitatively the evolution of 
the idealized pure-glue matter created in relativistic nu¬ 
clear collisions^. According to the QCD lattice calcula¬ 
tions (Celiketal. 1983a, 1983b; Karsch 2002; Borsanyi et 
al. 2012, Francis et al. 2015) this quarkless matter should 

This contradicts the conclusion of Gelis et al. (2004) that chemical 
nonequilibrium effects do not modify significantly the photon spectra. 

^ Such a case roughly coiTesponds to small rates of gg —» qq reactions. 
This can be simulated by choosing lai'ge r* within the approach developed 
in Sec. 2. However, we do not use now the ideal gas approximation for the 
gluonic matter. 



Fig. 6 Schematic picture of the temperature evolution of 
a high-energy collision in the pure glue scenario with the 
Yang-Mills first order phase transition to glueballs. 


T 



Fig. 7 (Color online) Transition temperature of baryon- 
free QGP versus the effective number of quark flavours. 


undergo the first order phase transition at the critical tem¬ 
perature Tc - 270 MeV. At this temperature the deconfined 
pure glue matter transforms into the confined state of the 
pure Yang-Mills theory, namely into a glueball fluid. 

Let us assume now that a hot thermalized gluon fluid 
(with no quarks and antiquarks) is created at the first stage 
of a relativistic collision. As the system cools and expands, 
it may reach the mixed phase domain at T — Tc and only 
after the glue plasma has completely transformed into the 
glueball fluid, the system cools down further. The system 
evolution in this pure SU(3) scenario is sketched in Fig.|3 
The possible appearance of super-cooled states and spinodal 
instabilities, associated with the first-order phase transition, 
can also be of particular interest. 

The heavy glueballs, produced in hadronization of 
a pure glue plasma, will evolve (presumably via cascade 
of two-body decays) into lighter states. Finally the system 
should decay into hadronic resonances and light hadrons. 
It was shown within the Frautschi - Hagedorn approach 
(Beitel et al. 2014) that the resulting yields of light hadrons 
and slopes of their spectra agree well with experimental data 
on heavy-ion collisions at RHIC and LHC energies. 
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In a more realistic scenario one should take into account 
that some quarks would already be produced before and dur¬ 
ing the first-order phase transition. Such a scenario could be 
modeled by introducing the time-dependent effective num¬ 
ber of quark degrees of freedom. With increasing number 
of quark degrees of freedom the temperature of the phase 
transition will decrease. At some point the first-order phase 
transition becomes a smooth crossover. This is schemati¬ 
cally shown in Fig. [T] The qualitative difference between 
the system evolution in the scenarios with the first-order 
and crossover transitions is demonstrated in Fig.|6] We think 
that the realization of a particular scenario depends on the 
energy and size of colliding objects. Thus, future studies of 
system-size and beam energy dependence of observables in 
nuclear collisions would be very useful. 

5 Conclusions 

The early stage of high multiplicity pp, pA and AA colli¬ 
sion events can represent a new state of deconfined mat¬ 
ter: a nearly quarkless, pure gluon plasma. According to the 
pure Yang - Mills lattice gauge theory, this matter under¬ 
goes, at a high temperature Tc - 270 MeV, the first-order 
phase transition into a confined Hagedorn-glueball fluid. 
Formation of such matter should lead to suppression of 
high pt photons and dileptons, to reduced baryon to meson 
ratios, as well as to enhanced yields of heavy (e.g. charmed) 
hadrons. We propose to search for signatures of pure glue 
states in the LHC/RHIC and cosmic rays experiments. 
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